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Abstract 

We re-examine the bounds on supersymmetric particle masses in light of the 
E821 data on the muon anomalous magnetic moment, a«. We confirm, extend and 
supersede previous bounds. In particular we find (at la) no lower limit on tan /3 
or upper limit on the chargino mass implied by the data at present, but at least 4 
sparticles must be lighter than 700 to 820 GeV and at least one sparticle must be 
lighter than 345 to 440 GeV. However, the E821 central value bounds tan/3 > 4.7 
and tuq < 690 GeV. For tan/3 < 10, the data indicates a high probability for 
direct discovery of SUSY at Run II or III of the Tevatron. 



Recently, the Brookhaven E821 experiment [JTJ] has reported evidence for a deviation 
of the muon magnetic moment from the Standard Model expectation. Immediately 
following that announcement appeared a number of papers analyzing the reported excess 
in terms of various forms of new physics, including supersymmetry (SUSY) [|2], || ||, ]5|, ||. 
The SUSY explanation is, for many of us, the most exciting of the various proposals 
since it implies SUSY at a mass scale not far above the weak scale. In particular, it 
implies a light slepton and a light gaugino, though "light" can still be as heavy as several 
hundred GeV. 

Previous analyses have generally concentrated on bounding the masses of the chargino 
and the sneutrino, or on the lightest of the sparticles. In this paper, we will re-examine 
the SUSY calculation, both in very general scenarios and in well-motivated, simplified 
SUSY models. We will show that there exist mass bounds on more than two SUSY 
particles, though which sparticles are bounded changes with tan/?. In particular, we 
will show that throughout parameter space, one may actually bound at least 4 spar- 
ticle masses if we are to use SUSY to explain the E821 data. The existence of these 
bounds will rely on very simple and clearly stated assumptions about the SUSY particle 
spectrum; these assumptions will not include a fine-tuning constraint. 



1 SUSY and 

The actual measurement performed by the E821 collaboration is of the muon's anoma- 
lous magnetic moment, which is to say, the coefficient a M of the non-renormalizable 
operator 



a 



2m M 

Within the Standard Model, a M is predicted to be |7j 

a^ M = 11 659 159.6(6.7) x 1(T 10 . 

It has been claimed that the quoted uncertainty is actually too small |J; the dominant 
uncertainty is the hadronic contribution to the photon polarization diagrams which is 
extracted from experimental measurements of R(e + e~ — > hadrons) in the vicinity of the 
low-energy meson resonances. However, strong arguments have been made to reinforce 
the quoted values ||, and we will accept them here as given. 
The measurement made by E821 is |J: 

ajf 21 = 11 659 202(14) (6) x 1(T 10 

from which one deduces a discrepancy between the experiment and the Standard Model 
of 

5a^ = 43(16) x 1(T 10 , 
that is, the measured value is larger than the prediction by 2.7a. 
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Figure 1: Super symmetric diagrams contributing to a M at one-loop. 



The SUSY contributions to a M have been known since the early days of SUSY and 



have become more complete with time[|H]]. In this paper we will follow the notation 
of Ref. [H which has the advantage of using the standard conventions of Haber and 



Kane |TT]; any convention which we do not define here can be found in either of these 
two papers. 

Before we begin the discussion of our work, let us review briefly a few of the analyses 
to date. The first paper to use the new E821 data in the context of SUSY was that of 
Czarnecki and Marciano who only attempted to approximately bound the sparticle 
spectrum. More complete analyses followed quickly thereafter by Everett et al. H and 
Feng and Matchev M. The former |3j argued that only very large values of tan/3 were 
generically consistent with the data; given those, a 1.5a upper bound of 450 GeV could 
be placed on some sparticle (gaugino or slepton) for tan (3 = 35, and 900 GeV for both 
a chargino and the muon sneutrino. The latter analysis Q found a model-independent 
la limit on the mass of the lightest "observable" sparticle of 490 GeV for tan (3 < 50. 
Martin and Wells did a more complete analysis, but focussed on the lighter chargino 
and charged smuon. Their model-independent analysis found no upper bound on the 
chargino mass at la but found < 500 GeV. With the added assumption of gaugino 
mass unification (i.e., M 2 — 2Mi), the lighter chargino was bounded by about 700 GeV at 
tan j3 = 30. The lack of a strong upper bound on the chargino was due to the neutralino 
contributions at small tan (3 (all the way down to tan (3 = 3) which can explain the 
E821 data without any chargino piece at all. Finally, there have been a number of other 
papers || which have studied the SUSY contributions within specialized scenarios which 
we will not discuss here. 



1.1 The diagrams 

In the mass eigenbasis, there are only two one-loop SUSY diagrams which contribute 
to a M , shown in Figure 1. The first has an internal loop of smuons and neutralinos, the 
second a loop of sneutrinos and charginos. But the charginos, neutralinos and even the 
smuons are themselves admixtures of various interaction eigenstates and we can better 
understand the physics involved by working in terms of these interaction diagrams, of 
which there are many more than two. We can easily separate the leading and subleading 
diagrams in the interaction eigenbasis by a few simple observations. 

First, the magnetic moment operator is a helicity-flipping interaction. Thus any 
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diagram which contributes to a M must involve a helicity flip somewhere along the fermion 
current. This automatically divides the diagrams into two classes: those with helicity 
flips on the external legs and those with flips on an internal line. For those in the first 
class, the amplitude must scale as m M ; for those in the second, the amplitudes can scale 
instead by mgusY, where msusY represents the mass of the internal SUSY fermion (a 
chargino or neutralino). Since wt-susy 3> Tn^, it is the latter class that will typically 
dominate the SUSY contribution to a M . Therefore we will restrict further discussion to 
this latter class of diagrams alone. 

Secondly, the interaction of the neutralinos and charginos with the (s)muons and 
sneutrinos occurs either through their higgsino or gaugino components. Thus each ver- 
tex implies a factor of either (the muon Yukawa coupling) or g (the weak and/or 
hypercharge gauge coupling). Given two vertices, the diagrams therefore scale as y 2 , 
gy^ or g 2 . In the Standard Model, y^ is smaller than g by roughly 10~ 3 . In the minmal 
SUSY standard model (MSSM) at low tan/?, this ratio is essentially unchanged, but 
because y^ scales as 1/ cos/3, at large tan/3 (~ 60) the ratio can be reduced to roughly 
10 _1 . Thus we can safely drop the y 2 contributions from our discussions, but at large 
tan (3 we must preserve the gy^ pieces as well as the g 2 pieces. [] 

The pieces that we will keep are therefore shown in Figures 2. In Fig. 2(a)-(e) are 
shown the five neutralino contributions which scale as g 2 or gy^; in Fig. 2(f) is the only 
chargino contribution, scaling as gy^. The contributions to from the ith neutralino 
and the mth smuon due to each of these component diagrams are found to be: 

g 1 N 2 1 X m iX m ,2 {BB) 

9ig2NiiN i2 X ml X m2 (WB) 

-V2g iyiI N a N i3 X 2 m2 (HB) 

-^giy.NMzX^ (BH) 

-^92y,N i2 N l3 X 2 ml (WH) 



1 m^m^ N 
A8n 2 m~. 



and for the /cth chargino and the sneutrino: 

1 777 777 " ~ ~ 

5a » = ~ 24tt 2 m 2 2< " Xk) 92V » Uk2Vkl - (WH) 



^ TV C 

The matrices N, U and V are defined in the appendix along with the functions F 2 ' . 
A careful comparison to the equations in the appendix will reveal that we have dropped 
a number of complex conjugations in the above expressions; it has been shown previ- 
ously |3], |I| that the SUSY contributions to a M are maximized for real entries in the mass 
matrices and so we will not retain phases in our discussion. 

In many of the previous analyses of the MSSM parameter space, it was found that it is 
the chargino-sneutrino diagram at large tan (3 that can most easily generate values of Sa^ 
large enough to explain the observed discrepancy. From this observation, one can obtain 
an upper mass bound on the lightest chargino and the muon sneutrino. However, this 

1 Pieces which are dropped from our discussion are still retained in the full numerical calculation. 
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Figure 2: Diagrams contributing to a M in the interaction eigenbasis. 

behavior is not completely generic. For example, Martin and Wells have emphasized that 
the BB neutralino contribution can by itself be large enough to generate the observed 
excess in a M , and since it has no intrinsic tan/3 dependence, they could explain the data 
with tan/3 as low as 3. We can reproduce their result in a simple way because the BB 
contribution has a calculable upper bound at which the smuons mix at 45°, <C m^, 
and Cm^^ with Ni = B. Then 



< - 38oo x io-- x ( m *> r ) ( 100GeV V 

I M | (BB) _ 327r2 ^ VIOOGeV/ V m Ai / 

where we have used the fact that (I11X12) < \ and < 3 and have included a 7% two- 
loop suppression factor. Though any real model will clearly suppress this contribution 
somewhat, this is still 10 2 times larger than needed experimentally. 

This pure BB scenario is actually an experimental worst-case, particularly for hadron 
colliders. The only sparticles that are required to be light are a single neutralino (which 
is probably 5-like) and a single /L The neutralino is difficult to produce, and if stable, 
impossible to detect directly. The neutralino could be indirectly observed in the decay of 
the jl as missing energy, but production of a /i at a hardon machine is highly suppressed. 
In the worst of all possible worlds, E821 could be explained by only these two light 
sparticles, with the rest of the SUSY spectrum hiding above a TeV. Further, even the 
"light" sparticles can be too heavy to produce at a 500 GeV linear collider. While this 
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case is in no way generic, it demonstrates that the E821 data by itself does not provide 
any sort of no-lose theorem for Run II of the Tevatron or even for the LHC and NLC. 

This raises an important experimental question: how many of the MSSM states must 
be "light" in order to explain the E821 data? In the worst-case, it would appear to be 
only two. Even in the more optimistic scenario in which the chargino diagram dominates 
5a M , the answer naively appears to be two: a single chargino and a single sneutrino. In 
this limit, 



\ Sa f\(CP) ^ 247r 2 



< 2600 x 10" 10 



x 



Ci 



VlOOGeVy V 



^100GeV\ 2 



m,- 



'tan p s 
30 



where we have bounded \F^\ by 10 by assuming ma < 1 TeV. But this discussion is 
overly simplistic and we can do much better, as we will see. 



1.2 Mass correlations 

There are a total of 9 separate sparticles which can enter the loops in Fig. H} 1 sneutrino, 
2 smuons, 2 charginos and 4 neutralinos. The mass spectrum of these 9 sparticles is 
determined entirely by 7 parameters in the MSSM: 2 soft slepton masses(mx, ttir), 2 
gaugino masses (Mi, M 2 ), the //-term, a soft trilinear slepton coupling (Ap) and finally 
tan j3. Of these, Ap plays almost no role at all and so we leave it out of our discussions 
(see the appendix). And in some well- motivated SUSY-breaking scenarios, Mi and M 2 
are also correlated. Thus there are either 5 or 6 parameters responsible for setting 9 
sparticle masses. There are clearly non-trivial correlations among the masses which can 
be exploited in setting mass limits on the sparticles. 

First, there are well-known correlations between the chargino and neutralino masses; 
for example, a light charged Cj ~ W implies a light neutral Nj ~ W and vice-versa. 

There are also correlations in mixed systems (i.e., the neutralinos, charginos and 
smuons) between the masses of the eigenstates and the size of their mixings. Consider 
the case of the smuons in particular; their mass matrix is given in the appendix. On 
diagonalizing, the left-right smuon mixing angle is given simply by: 

2m u utan/3 
tan 29 n ~ ^ 



11 ~ Ml -M^' 

The chargino contribution is maximized for large smuon mixing and large mixing occurs 
when the numerator is of order or greater than the denominator; since the former is 
suppressed by m M , one must compensate by having either a very large //-term in the 
numerator or nearly equal Ml and Mr in the denominator, both of which have profound 
impacts on the spectrum. 

There is one more correlation/constraint that we feel is natural to impose on the 
MSSM spectrum: slepton mass universality. It is well-known that the most general 
version of the MSSM produces huge flavor-changing neutral currents (FCNCs) unless 
some external order is placed on the MSSM spectrum. By far, the simplest such order is 
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the demand that sparticles with the same gauge quantum numbers be degenerate. Thus 
we expect m fL = = thl and rrif R = m^ R = thr. Then the mass matrix for the 
stau sector is identical to that of the smuons with the replacement m M — > m T in the off- 
diagonal elements. This enhancement of the mixing in the stau sector by m T /m^ ~ 17 
implies that m fl < . In particular, if 



then < and QED will be broken by a stau vev. Given slepton universality, this 
imposes a constraint on the smuon mass matrix: 



where Ml ^ are the positive roots of M| R . While not eliminating the possibility of 
9^ ~ 45°, this formula shows that a fine-tuning of at least 1 part in 17 is needed to 
obtain 0(1) mixing. We will not apply any kind of fine-tuning criterion to our analysis, 
yet we will find that this slepton mass universality constraint sharply reduces the upper 
bounds on slepton masses which we are able to find in our study of points in MSSM 
parameter space. 

(As an aside, if one assumes slepton mass universality at some SUSY-breaking mes- 
senger scale above the weak scale, Yukawa-induced corrections will break universality by 
driving the stau masses down. This effect would further tighten our bounds on smuon 
masses and mixings.) 

The above discussion has an especially large impact on the worst-case scenario in 
which the BB contibutions dominates 5a^. For generic points in MSSM parameter 
space, one expects that tan 26 < 1/17 which reduces the size of the BB contribution by 
a factor of 17. As a byproduct, the masses required for explaining the E821 anomaly 
are pushed back towards the range that can be studied by a 500 GeV linear collider. 

2 Numerical results 

Now that we have established the basic principle of our analysis, let us carry it out in 
detail. We will concentrate on three basic cases. The first case is the one most often 
considered in the literature: gaugino mass unification. Here one assumes that the weak- 
scale gaugino mass parameters (Mi and M 2 ) are equal at the same scale at which the 
gauge couplings unify. This implies that at the weak scale M x = (5/3)(a 1 /a 2 )M 2 . The 
second case we consider is identical to the first with the added requirement that the 



M 2 L M 2 R < m 2 /i 2 tan 2 (3 




or on the smuon mixing angle: 
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lightest SUSY sparticle (LSP) be a neutralino. This requirement is motivated by the 
desire to explain astrophysical dark matter by a stable LSP. Finally we will also consider 
the most general case in which all relevant SUSY parameters are left free independent 
of each other; we will refer to this as the "general MSSM" case. 

The basic methodology is simple: we put down a logarithmic grid on the space of 
MSSM parameters (Mi, M2, m^, ttlr and fi) for several choices of tan/?. The grid 
extends from 10 GeV for Mi, M 2 and fi, and from 50 GeV for m^, up to 2 TeV for all 
mass parameters. For the case in which gaugino unification is imposed M 2 is no longer 
a free parameter and our grid contains 10 8 points. For the general MSSM case our grid 
contains 3 x 10 9 points. Only fi > is considered since that maximizes the value of 5a^. 
Finally, for our limits on tan/3 we used an adaptive mesh routine which did a better job 
of maximizing over the space of MSSM inputs. By running with grids of varying 
resolutions and offsets we estimate the error on our mass bounds to be less than ±5%. 

2.1 Bounds on the lightest sparticles 

Perhaps the most important information that can be garnered from the E821 data is 
an upper bound on the scale of sparticle masses. In particular, one can place upper 
bounds on the masses of the lightest sparticle(s) as a function of 5a^. Previous analyses 
have often followed this approach, deriving upper bounds on the lightest sparticle from 
among the gauginos and sleptons, or even more specifically, from among the charginos 
and smuons. We too derive bounds on the lightest sparticles, but as we have argued in 
the previous section, we can also derive bounds on the second, third, and even fourth 
lightest sparticles through the mass correlations. 

These bounds on additional light sparticles provide an important lesson. Without 
them there remains the very real possibility that the E821 data is explained by a pair of 
light sparticles and that the remaining SUSY spectrum is out of reach experimentally. 
But our additional bounds will give us some indication not only of where we can find 
SUSY, but also of how much information we might be able to extract about the fun- 
damental parameters of SUSY — the more sparticles we detect and measure, the more 
information we will have for disentangling the soft-breaking sector of the MSSM. 

In Fig. ^| we have shown the upper mass bounds for the lightest four sparticles 
assuming gaugino mass unification. These bounds are not bounds on individual species 
of sparticles (which will come in the next section and always be larger than these bounds) 
but simply bounds on whatever sparticle happens to be lightest. The important points 
to note are: (i) the maximum values of the mass correspond to the largest value of 
tan p, which is to be expected given dominance of the chargino diagram at large tan /3; 
(ii) the la limit (central value) of the E821 data requires at least 4 sparticles to lie 
below roughly 700 (500) GeV; and (in) for low values of tan/3 a maximum value of Sa^ 
is reached (we will return to this later). 

The same plots could be produced with the additional assumption that the LSP 
be a neutralino, but we will only show the case for the LSP bound, in Fig. In this 
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Figure 3: Bounds on the masses of the four lightest sparticles as a function of Sa^ for tan f3 = 3, 
5, 10, 30 and 50. These figures assume gaugino unification only. 
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Figure 4: Bound on the mass of the LSP as a function of Sa^ for tan/3 = 3, 5, 10, 30 and 
50. The dotted lines assume gaugino unification only while the solid lines require additionally 
that the LSP is a neutralino. The dashed line is the bound in the general MSSM, calculated 
at tan/3 = 50. 



figure, the solid lines correspond to a neutralino LSP, while the dotted lines are for 
the more general case discussed above (i.e., they match the lines in Fig. 0(a)). Notice 
that for Sa^ > 40 x 10 -10 there is little difference between the cases with and without 
a neutralino LSP. Furthermore, at the extreme upper and lower values of tan/3 there 
is little difference. It is only for the intermediate values of tan/3 that the mass bound 
shifts appreciably; for tan /3 = 10 it comes down by as much as 50 GeV compared to the 
more general case. 

Finally, we consider the most general MSSM case, i.e., without gaugino unification. 
Here the correlations are much less pronounced, but interesting bounds still exist. For 
example the central value of the E821 data still demands at least 3 sparticles below 
500 GeV (rather than four for the previous cases). In figure |5| we demonstrate this 
explicitly by plotting the masses of the four lightest sparticles for tan/3 = 50 and a wide 
range of 5a^. We see that dropping the gaugino unification requirement has one primary 
effect: the mass of the LSP is significantly increased. This is because the LSP in the 
unified case is usually a Aq ~ B but isn't itself responsible for generating Sa^. In the 
general case, the LSP must participate in Sa^ (otherwise its mass could be arbitrarily 
large) and so is roughly the mass of the second lightest sparticle in the unified case, 
whether that be a jl or C. Otherwise the differences between the more general MSSM 
and the gaugino unified MSSM are small. In particular we still find that at least 4 
sparticles must be light, though the lex bound of 700 GeV for the unified case extends 
now slightly to 820 GeV. 
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Figure 5: Bounds on the masses of the four lightest sparticles as a function of 5a „ for tan (5 = 50. 
The dotted lines assume gaugino unification only while the solid lines are for the general MSSM. 

We have summarized all this data on the LSP in Table [I] where we have shown the 
mass bounds (using both the la limit and the central value of the E821 data) on the 
LSP for various tan/5 values and with our various assumptions. The last line in the 
table represents an upper bound for any model with tan/3 < 50: tt^lsp < 440 (345) GeV 
for the E821 la lower bound (central value) of Sa^. But perhaps of equal importance 
are the bounds on the next 3 lightest sparticles (the "2LSP," "3LSP" and "4LSP"): 
m 2LSP < 460 (355) GeV, m 3LSP < 600 (465) GeV and m 4LS p < 820 (580) GeV. Thus 
there must be at least 4 sparticles below 820 GeV even in the most general MSSM 
scenario, and at least two sparticles accessible to a \/s = 1 TeV linear collider. 

2.2 Bounds on the sparticle species 

In the previous subsection, we derived bounds on the lightest sparticles, independent 
of the identity of those sparticles. Another important piece of information that can be 
provided by this analysis is bounds on individual species of sparticles, for example, on 
the charginos or on the smuons. These bounds will of mathematical necessity be higher 
than those derived in the previous section, but still provide important information about 
how and where to look for SUSY. In particular, they can help us gauge the likelihood 
of finding SUSY at Run II of the Tevatron or at the LHC. 

There is one complication in obtaining these bounds. At low tan f3 the data is most 
easily explained by the neutralino diagrams and as such there must be at least one light 
smuon and one light neutralino. At larger tan/5 values of Sa^ as large as the E821 
central value generally require contributions from the chargino diagrams, so there must 
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440 (345) 




345 f265j 
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Table 1: Upper bounds on the mass (in GeV) of the lightest sparticle for the general MSSM, 
the MSSM assuming gaugino mass unification, and the MSSM with gaugino mass unification 
plus a neutralino LSP. The entries represent the bound for the la limit (central value) of the 
E821 data. The boldfaced tan/3 = 50 entries represent upper bounds over all tan (5 < 50. 
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Figure 6: Bounds on the masses of N\ and /2i as a function of 5a^ for various tan (3 with 
gaugino unification assumed. 

be a light chargino and a light sneutrino. However the correlations already discussed 
preserve the bounds on the various species over the whole range of tan (3. A bound on 
frij) implies a bound on m^, and a bound on mg implies a bound on at least one of 
the rrifj,, and in certain cases (such as gaugino unification), the converses may be true 
as well. 

We have shown in Fig. || the mass bounds on fti and N% under the assumption of 
gaugino unification; a plot for C1/N2 will appear later in our discussion of Tevatron 
physics. Note that N\ must be relatively light, even for large tan/?, thanks to the 
gaugino unification condition, while jli can be heavier but must still lie below 820 GeV 
at lcr. 

Finally, we can consider the general MSSM without gaugino unification. The results 
can best be summarized by showing the mass bounds on the various sparticles at tan j3 = 
50 in relations to their bounds in the unified case. This is done in Fig. [?]. We can see 
from the figure that the bound on the fli is essentially identical to that in the gaugino 
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Figure 7: Bounds on the masses of fix, C\, and N\ as a function of 5a^ for tan/3 = 50, 
with gaugino unification (solid) and in the general MSSM (dashed). The two lines for the /ii 
essentially overlap. 

unification picture. However the gaugino masses have shifted, and the reason is no 
mystery. Once again, as discussed in the previous section, the lightest neutralino is 
no longer a B-\ike spectator to the magnetic moment, but is a VF-like partner of a 
participating VF-like chargino. 

We summarize our results for the various sparticle species in Table |2|. There we 
have shown the upper bounds on several sparticles in the general MSSM, the MSSM 
with gaugino unification, and the previous case with the additional requirement of a 
neutralino LSP ("dark matter"). The bounds represents those obtained using the lex 
limit (central value) of the E821 data. 

2.3 Bounds on tan j3 

The final bound we will derive using the E821 data is on tan/3. There has been some 
discussion in the literature about which values of tan /3 are capable of explaining the 
data. And in fact, we concur that at lower tan/3, there is a real suppression in the 
maximum size of 5a fl . So in Figure ^| we have shown the maximum attainable value 
of Sa^ as a function of tan /3 with and without the added assumption of gaugino mass 
unification. 

The limit in Fig. |8| clearly divides into two regions. At 5a^ > 36 x 10~ 10 the chargino 
contribution dominates and thus 5a M oc y^, scaling linearly with tan/3. At lower Sa^, 
however, both neutralino and chargino contributions can be important so it becomes 
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470 (350) 
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830 (565) 


575 (440) 
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830 (565) 
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870(680) 


825 (665) 


825 (665) 
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865 (675) 


820 (660) 
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Table 2: Upper bounds on the mass (in GeV) of various sparticles for the general MSSM, the 
MSSM assuming gaugino mass unification, and the MSSM with gaugino mass unification plus 
a neutralino LSP. The entries represent the bound for the la limit (central value) of the E821 
data. These bounds are for all tan/3 < 50. 
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Figure 8: Bounds on tan (3 as a function of <5a M for the general MSSM. The dotted lines 
represent the E821 central value and ±lcr bounds on da^. 
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possible to generate 5a^ with much smaller values of tan (3 than would be possible from 
the charginos alone. The central value of the E821 data implies tan/? > 4.7; but already 
at la all values of tan f3 down to 1 are allowed. This result contradicts most statements 
made in the literature to date about limiting tan/5 using a^. at present we are not 
able to place any la bound on tan/3 larger than 1 due to the very real possibility of 
a neutralino-dominated 5a M . Future reduction in the error bars on a M will provide an 
important opportunity for placing a meaningful lower bound on tan /3. 



2.4 Implications for the Tevatron 

At its simplest level, the measurement of £a M , an anomaly in the lepton sector, has little 
impact on the Tevatron, a hadron machine. In particular, the light smuons associated 
with 5a^ cannot be directly produced at the Tevatron, occuring only if heavier (non- 
leptonic) states are produced which then decay to sleptons. In the calculation of a M , 
the only such sparticles are the neutralinos and charginos. These states can be copi- 
ously produced and in fact form the initial state for the "gold-plated" SUSY trilepton 
signature. 

Of particular interest for the trilepton signature are the masses of the lighter chargino 
(Ci) and 2nd lightest neutralino (N 2 ). Studies of mSUGRA parameter space indicate 
that the sensitivity to the trilepton signature at Run II/III of the Tevatron depends 
strongly on the mass of sleptons which can appear in the gaugino decay chains. For heavy 
sleptons, the Tevatron is only sensitive to gaugino masses in the range JT2] tuq ^ < 130 



to 140 GeV for 10 ffcT 1 of luminosity and 145 to 155 GeV for 30 fb _1 , with smaller values 
for smaller tan (3. However, for light sleptons (below about 200 GeV) the range is con- 
siderably extended, up to gaugino masses around 190 to 210 GeV. Of course, this latter 
range is exactly the one most relevant for understanding a^. 

It is impossible in the kind of analysis presented here to comment on the expected 
cross-sections for the neutralino-chargino production (there is no information in a M on 
the masses of the t-channel squarks, for example) but we can examine the mass bounds 
on C\ and N 2 . In Fig. || we have shown just that: the upper bound on the heavier of 
either C\ or N 2 as a function of Sa^ for several values of tan j3. 

A few comments are in order on the plot. First, the plot assumes gaugino unification; 
dropping that assumption can lead to significantly heavier masses for the N 2 though not 
for the Ci, this is because the light C\ is needed to participate in the magnetic moment, 
while N 2 can decouple. Second, we have also assumed a neutralino LSP; this is to be 
expected since the event topology for the trilepton signal assumes a stable, neutralino 
LSP. Finally, on the y-axis is actually plotted , but in every case we examined with 
gaugino unification, the difference in the maximum masses of C\ and N 2 differed by at 
most a few GeV. This is because they are both dominantly wino-like in the unified case 
and thus have masses ~ M 2 . 

From the figure it is clear that a no-lose theorem for the Tevatron is not lurking in 
the current E821 data. However, if the central value reported by E821 holds up and 
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Figure 9: Mass bounds on C\ and N% (where — rrifj ) as a function of 5a^ for tan/? = 3, 
5, 10, 30 and 50. The dotted lines represent the E821 central value and la bounds on <5a M . 
This figure assumes gaugino unification and a neutralino LSP. 



tan (3 < 10 then one should expect the Tevatron to find a trilepton signal for SUSY. So 
there is actually significant hope for a positive result. We cannot emphasize enough too 
that these are upper bounds on the sparticle masses and in no way represent best fits or 
preferred values. Thus even for larger tan/3 or smaller Sa^, there is good reason to hope 
that the Tevatron will be able to probe the gaugino sector in Run II or III. 



2.5 Implications for a Linear Collider 

A concensus is currently forming in favor of building a y/s = 500 GeV linear collider, 
presumably a factory for sparticles with masses below 250 GeV. What does the mea- 
surement of a M tell us with regards to our chances for seeing SUSY at y/s = 500 GeV? 
And how many sparticles will be actually accessible to such a collider? 

The analysis of the previous section can put a lower bound on the number of ob- 
servable sparticles at a linear collider as a function of Sa^ and tan (3 and we show those 



numbers as a histogram in Fig. [It]. In this figure, we have shown the minimum number 
of sparticles with mass below 250 GeV for tan/3 = 5, 10, 30 and 50, assuming gaugino 
unification. In the graph, the thinner bars represent smaller tan {3. As is to be expected, 
the number of light states increases with increasing 5a^ and with decreasing tan (3. How- 
ever note that there are no tan/3 = 5 lines for Sa^ > 40 x 10~ 10 since there is no way to 
explain such large <5a M values at low tan (3. 

We see from the figure also that for tan (3 > 30 there is no guarantee that a 500 GeV 
machine would produce on-shell sparticles; this is not to be taken to mean that one 
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Figure 10: Minimum number of sparticles directly observable at a y/s = 500 GeV linear collider 
as a function of Sa^. The bars represent tan (3 = 5, 10, 30 and 50 in order of increasing thickness. 
This graph assumes gaugino unification but does not include additional sleptons implied by 
slepton mass universality. 

should not expect their production, simply that cannot guarantee it. However for 
tan ft < 10, the la limit on a M from E821 indicates that at least 2 to 3 sparticles will 
be accessible to a 500 GeV machine. This counting does not include extra sleptons due 
to slepton mass universality; for example, a light muon sneutrino also implies light tau 
and electron sneutrinos, and likewise for the charged smuon. 

A similar bar graph can also be made for a 1 TeV machine, though we do not show 
it here. However the relevant numbers can be inferred from Fig. 0; we see that such a 
machine is guaranteed to produce at least three sparticles for tan ft < 50 at the E821 la 
limit. 

3 Conclusions 

Deviations in the muon anomalous magnetic moment have long been advertised as a key 
hunting ground for indirect signatures of SUSY. And it may be that we finally have the 
evidence we need in the 2.7a deviation reported by E821. If this signal is real and if 
SUSY is the correct explanation for it, then light sparticles are requisite. In this analysis, 
we have confirmed, extended and overridden some of the bounds on these light sparticles 
that have appeared previously. In particular, for the la lower bound on 5a^ derived from 
the E821 data, we obtain the following for the most general MSSM parametrization: 

• there must be at least 4 sparticles with masses below 820 GeV (700 GeV for unified 
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gauginos) ; 

• the lightest sparticle must lie below 440 GeV (345 GeV for unified gauginos); 

• there is no lower bound on tan/5 (the central value of the E821 data requires 
tan/3 > 4.7); 

• the fx\ and the v must lie below 870 GeV while one neutralino must fall below 
610 GeV (those bounds become 680 and 455 GeV respectively for unified gauginos); 

• there is no upper bound on mg (the data's central value requires < 690 GeV); 

• if tan (3 < 10, then a 500 GeV linear collider is guaranteed to produce at least 2 
sparticles, and Run II /III of the Tevatron is likely to see a trilepton signature 
(assuming a stable neutralino LSP); for tan/3 > 10, neither of these can be guar- 
anteed. 

In the most interesting theoretical scenario, i.e. gaugino unification with a neutralino 
LSP, the bounds on SUSY masses can be much lower, with the neutralino LSP below 
340 GeV, C\ and N2 below 575 GeV and sleptons below 825 GeV. Our only important 
simplifications were to drop the dependence on Ap_ from the smuon mixing matrices, and 
to assume slepton mass universality. 

If E821 has discovered evidence for SUSY, then it is clear that we have many more 
discoveries awaiting us. The current data implies that the LHC will find SUSY directly; 
such a discovery for the Tevatron or the NLC is not guaranteed, though the data on 
greatly enhances the likelihood that these programs will be successful. 
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Appendix 

The supersymmetric contributions to a M are generated by diagrams involving charged 
smuons with neutralinos, and sneutrinos with charginos. The most general form of the 
calculation, including phases, takes the form (we follow Ref. ||): 

i,m \ ftm Am / 
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Sap 



TO, 



16tt 2 ^ \ 12m; 



m "j(l4l 2 + l<£ | 2 )if + ^Re[cjcf]if (*») 



3m~ 



(1) 



where % = 1,2,3,4, to = 1,2, and fc = 1,2 label the neutralino, smuon and chargino 
mass eigenstates respectively, and 

nf m = VlgiNaX^ + y„N i3 X ml , 



= —(g 2 N i2 +g 1 N il )X* ml -y^N i3 X* m2 , 

= y^,u k2 , 



(2) 



Vii = fi^^v/ V2mw cos (3 is the muon Yukawa coupling, and g\ >2 are the U(l) hypercharge 
and SU(2) gauge couplings. The loop functions depend on the variables x im = m^Jmj^, 
Xk = rrig / to- and are given by 
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F 2 c (x) 
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1 - Qx + 3x 2 + 2x 3 - Qx 2 lux 

1 — x 2 + 2x In x 
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x) 



(3) 



F 2 N (1) 



For degenerate sparticles (x = 1) the functions are normalized so that F^(l) 
F^\l) = F 2 (l) = 1. We can also bound the magnitude of some of these functions; in 
particular \F 2 N (x)\ < 3 while |F 2 c (x)| is unbounded as x — > 0. 
The neutralino and chargino mass matrices are given by 

/ 



Mi 
M 2 

mzswcp mzcwcp 

\ m z s w s /3 -m z c w sp 



-mzswcp m z s w Sj3 \ 
m z cwC/3 -mzcwsp 





-fj, 



-n 





(4) 



/ 



and 



fx 



(5) 



M 2 
^ y/2m w Cp 

where sp = sin/5, cp = cos/3 and likewise for 8w- The neutralino mixing matrix Nij and 
the chargino mixing matrices U k i and V k i satisfy 



N*MfjN ] = diag^^TO^TO^TO^J 
U*M d V* = diag(m^,m^). 



(6) 
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The smuon mass matrix is given in the {ft,L, P>r} basis as: 

M 2_( M 2 L m^At-fMttmP) \ 

where 

Ml = m 2 L + (s 2 w -^)m 2 z cos2f3 

Ml = m| - s 2 w m 2 z cos 2(3 (8) 

for soft masses m 2 L and m\\ the unitary smuon mixing matrix X mn is defined by 

XM*X* = diag (ml, ml). (9) 

We will define a smuon mixing angle 9p, such that Xu = cos9p, and X12 = sin6^. In our 
numerical calculations we will set = 0. At low tan (3 we have checked that varying 
Ap makes only slight numerical difference, while at large tan/3 it has no observable effect 
whatsoever. Finally, the muon sneutrino mass is related to the left-handed smuon mass 
parameter by 

m~ = m 2 L + cos 2/5- (10) 



The leading 2-loop contributions to Sa^ have been calculated |13| and have been 
found to suppress the SUSY contribution by a factor (ia/ir) log(m SU sY l m l) ~ 0.07; we 
will include this 7% suppression in all our numerical results. 
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